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A method and apparatus which allow accurate spectrophotometric determinations of the concentrations of various hemog- 
lobin species in whole blood without hemolysis or dilution. To overcome complex optical properties of whole blood, the invention 
employs 1) an optical apparatus (10, 1 1, 12, 13 and 14), designed to maximize the true optical absorbance of whole blood and to 
minimize the effects of tight scattering on the spectrophotometric measurements of concentrations of various constituent compo- 
nents, and 2) methods to correct the hemoglobin concentration measurements for light scattering and for the effects of the finite 
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cal parameters, such as sample thickness, detector size and shape, sample to detector distance, wavelengths, monochromicity, 
maximum angle of light capture by detector, are optimal values to minimiz e the contribution of light scattering to the total optical 
attenuation of unaltered whole blood and maximize contribution of true optical absorbance. 
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DESCRIPTION 

METHOD AND APPARATUS FOR DIRECT SPECTROPHOTOMETRY 
MEASUREMENTS IN UNALTERED WHOLE BLOOD 

5 A portion of the disclosure of this patent document 

contains material which is subject to copyright 
protection. The copyright owner has no objection to the 
facsimile reproduction by anyone of the patent document 
or the patent disclosure, as it appears in the Patent and 
10 Trademark Office patent file or records, but otherwise 
reserves all copyright rights whatsoever. 



FIELD OF THE INVENTION 



15 This invention relates to a method and apparatus to 

assess the optical transmittance of a sample of unaltered 
whole blood at multiple wavelengths to attain an accurate 
measurement of its total hemoglobin concentration, the 
concentration of bilirubin, and the concentrations of 

20 oxy-, deoxy- , carboxy-, met-, and sulf hemoglobin. 



To the best of the Applicants' knowledge, no prior 
25 art whether patented or not has ever successfully 

exploited the optical transmittance of unaltered whole 
blood to achieve an accurate measurement of the total 
hemoglobin concentration (THb) in a blood sample that 
could possibly contain as many as five different species 
30 of hemoglobin. Several different properties of whole 
blood have precluded such measurements. These 
characteristics include 1) the complex chemical mixture 
of as many as five individual hemoglobin species that are 
present in varying concentrations in a typical blood 
35 sample, and 2) the complex optical properties of whole 
blood that result from a combination of true optical 
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absorbance by multiple light-absorbing compounds 
(including the hemoglobins , bilirubin, and other 
pigments) , and the pronounced light-scattering ability of 
the red blood cells that are present in high 
5 concentrations. 

Many prior methods have been devised to measure the 
total hemoglobin concentration in the complex chemical 
mixture of oxy-, carboxy- , met-, deoxy-, and 

10 sulf hemoglobin that comprise the total hemoglobin 

composition. For example, Loretz (U.S. Patent 4,357,105) 
teaches a method whereby the hemoglobin is released from 
the red blood cells by the process of hemolysis, i.e. by 
disrupting the red blood cell membrane. Simultaneously, 

15 the various hemoglobin species are diluted by reagent 

solutions that convert them into a single chemical form, 
cyanmethemoglobin, the concentration of which is then 
measured spectrophotometrically by conventional 
apparatus . 

20 

Figure 1 presents a graph of the extinction 
coefficients (indicative of optical absorbance) as a 
function of radiation wavelength for five individual 
hemoglobin species. Because the individual hemoglobin 

25 species differ significantly from each other in their 
optical absorbance spectra, no method that employs a 
single wavelength can measure the total hemoglobin 
concentration spectrophotometrically. Therefore, the 
chief advantage of the many prior methods that convert 

30 the sundry species into a single known species is that 
only one wavelength is needed for the subsequent 
spectrophotometric measurement of the resulting single 
species. The disadvantages of the above-mentioned 
chemical conversion methods are that 1) hemolysis and 

35 accurate dilutions are needed, 2) the reagents are often 
toxic, 3) chemical conversion is slow, 4) the red blood 
cells and the chemical nature of the whole blood sample 
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are so drastically altered that the same sample cannot be 
used for further subsequent analyses of other blood 
constituents, and 5) the chemical conversion of the 
various hemoglobin derivatives into a single species 
5 makes it impossible to measure the original 

* 

concentrations of the several individual species. The 
latter disadvantage is a serious limitation because the 
concentrations of oxy-, carboxy-, met-, deoxy-, and 
sulfhemoglobin provide valuable diagnostic information in 
10 many different medical applications. 

In several patented methods (U.S. Patents 4,134,678; 
4,013,417; 3,972,614; and 4,997,769), the need for 
chemical conversion of the various hemoglobin derivatives 

15 into a single hemoglobin species has been eliminated by 
using multiple appropriate wavelengths to measure the 
concentration of each species. In such methods, the 
blood sample is either hemolyzed and diluted (Brown et 
al., U.S. Patent 4,134,678; Raffaele, U.S. Patent 

20 4,013,417) or just hemolyzed (Johansen et al., U.S. 

Patent 3,972,614; Lundsgaard, U.S. Patent 4,997,769) to 
eliminate the light scattering by the red blood cells. 
The hemolyzed sample is then irradiated with at least one 
distinct wavelength for each of the three, four, or five 

25 hemoglobin species presumed to be present in the sample. 
By measuring the optical density of the hemolyzed sample 
with at least as many different wavelengths as there are 
unknown individual hemoglobin species, such methods yield 
a set of simultaneous equations that can subsequently be 

30 solved for the relative concentration of each of the 
various hemoglobin derivatives. These prior methods 
eliminate the need for chemically converting the various 
hemoglobin derivatives into a single species, but they 
suffer from the disadvantage of requiring a complex, 

35 bulky, and expensive apparatus comprised of pumps, 
plumbing, associated control circuitry, and in some 
cases, ultrasonic hemolyzers to dilute and hemolyze the 
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blood sample. Finally, the techniques that employ 
hemolysis have two additional disadvantages: l) their 
complicated plumbing systems are prone to clogging by 
blood residue, and 2) they aspirate and destroy the blood 
5 sample so that it cannot be retrieved and subjected to 
further analysis of other constituents. 

In whole blood, the major obstacle to making optical 
measurements is the intense light scattering caused by 

10 the highly concentrated red blood cells, e.g. 5.4 X 10 6 
RBCs//il for human males. Consequently, all prior art has 
relied exclusively on hemolysis or on both hemolysis and 
dilution to reduce or eliminate the light scattering 
caused by red blood cells. This light scattering arises 

15 because the blood plasma and the red blood cells have 

different indices of refraction (B. Bearer and S. Joseph, 
Quart. J. Hicrosc. Sci. 95:399,1954). Many 
uncontrollable factors make such light scattering 
unpredictable from one blood sample to another, including 

20 1) the different plasma protein concentrations that 
determine the refractive index of plasma, 2) the 
aggregation and sedimentation of red blood cells in the 
sample, 3) the different hemoglobin concentrations inside 
the red blood cells that alter their refractive index, 

25 and 4) the size and shape of the red blood cells. The 
size and shape of the red blood cells differ from one 
animal species to another; they are also altered by 
disease states (e.g. sickle-cell anemia and 
spherocytosis) ; and they are altered by the osmotic 

30 pressure of the blood. Thus all prior methods which 
measure the total hemoglobin concentration and the 
relative concentrations of individual hemoglobin species 
have employed hemolysis to eliminate this intense, 
unpredictable light scattering caused by the red blood 

35 cells. 
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Prior art that relies on hemolysis to eliminate the 
intense light scattering by the red blood cells has the 
further disadvantage that even after thorough hemolysis 
the sample can be relatively turbid (L.R. Sehgal et al., 
5 Critical Care Medicine, 12:907-909, 1984). This residual 
turbidity in hemolyzed blood is small compared with the 
intense light scattering of red blood cells in unaltered 
whole blood (J.M. Steinke and A. P. Shepherd, IEEE Trans . 
Biomed . Eng . BME-33: 294-301, 1986). Nevertheless, the 

10 residual turbidity in hemolyzed blood causes troublesome 
errors in the hemoglobin measurements. Thus, an 
additional disadvantage of prior apparatus that rely on 
hemolysis is the residual turbidity that results from a 
small number of unlysed red blood cells, lipid particles 

15 such as chylomicrons (a normal constituent of plasma that 
persists after hemolysis) , light- scattering cell 
fragments produced by the hemolysis process, and other 
causes that are unknown. Although Lundsgaard (U.S. 
Patent 4,997,769) teaches a method for dealing with the 

20 residual turbidity that remains after the hemolysis 

process, it would be advantageous to eliminate hemolysis 
completely and make spectrophotometric measurements 
directly in unaltered whole blood. 



25 The Applicants have found that the scattering of 

light by unaltered (unhemolyzed) whole blood differs in 
five crucial ways from the turbidity of hemolyzed blood. 
First, as mentioned previously, the turbidity of 
hemolyzed blood is insignificant in magnitude in 

30 comparison with the light scattering of unhemolyzed whole 
blood. This conclusion is reinforced by Figure 2 which 
shows a measure of light scattering in whole blood before 
and after hemolysis. For these data, the light 
scattering in unaltered, whole blood is approximately 20 

35 times greater than after hemolysis. In fact, light 

scattering by red blood cells in unaltered whole blood is 
so intense that scattered light predominates over 
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unscattered light after it has passed through as few as 
10 layers of red blood cells (C.C. Johnson. J. Assoc. 
Adv. Med. Instrument. 4:22-27,1970). Prior art was not 
designed to accommodate the greater magnitude of the 
5 scattering effects of unaltered whole blood. 



Second, a fundamental observation is that red blood 
cells scatter light at relatively large angles (J.M. 
Steinke and A. P. Shepherd, Applied Optics 27:4027- 

10 4033,1988). The Applicants have discovered that, even 
though most of the light is scattered at small angles, 
the magnitude of large-angle light scattering is 
sufficient to cause serious errors in spectrophotometric 
measurements on whole blood. Prior art does not address 

15 the problem of designing practical instruments by 

capturing the large-angle light scattering of unhemolyzed 
blood. 



Third, the light scattering by red blood cells 
20 varies with wavelength in a complicated manner due to 

physical factors such as the relative size of red blood 
cells with respect to the wavelength, and the shape of 
the red blood cells. Figure 3 shows the wavelength 
dependence of light scattering by RBCs as predicted by 
25 optical theory (R.N. Pittman and B.R. Duling, J. Appl. 

Physiol. 38:315-320,1975; R.N. Pittman, Ann. Biomed. Eng. 
14:119-137,1986). Prior art does not address the complex 
wavelength dependence of the scattering effects of 
unaltered whole blood. 

30 

Fourth, as Figure 4 shows, the light scattering by 
red blood cells is not only a function of wavelength, but 
it also depends on the composition of the hemoglobin 
contained in the red blood cells, i.e. the light 
35 scattering depends on the relative concentrations of 

oxy-, deoxy-, carboxy- , met-, and sulfhemoglobin in the 
red blood cells. Prior art does not address the complex 
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dependence of light scattering on the concentrations of 
the individual hemoglobin species in unaltered whole 
blood . 

5 Fifth, many poorly understood, uncontrolled 

processes occur in unaltered, whole blood that change its 
optical properties. Such miscellaneous causes of 
unpredictable light scattering include 1) the different 
plasma protein concentrations that determine the 

10 refractive index of plasma in one sample versus another, 
2) the aggregation of red blood cells in the sample, 3) 
the different hemoglobin concentrations inside the red 
blood cells that alter their refractive index, 4) the 
size and shape of the red blood cells, 5) chylomicrons or 

15 other light-scattering lipid particles, 6) cell 

fragments, 7) microscopic clots, and 8) light-sieving 
effects of sedimented RBCs. 

As shown in the comparative example below, the 
20 apparatus of Lundsgaard, U.S. patent 4,997,769, fails to 
yield valid results when measuring a sample of unaltered 
(unhemolyzed) whole blood. Similarly, other prior art 
that relies on hemolysis such as Brown et al., U.S. 
Patent 4,134,678; Raffaele, U.S. Patent 4,013,417; and 
25 Johansen et al., U.S. Patent 3,972,614 would also fail to 
yield valid results on unhemolyzed whole blood because 
these methods also fail to capture the light scattered at 
large angles by the sample, and because their 
measurements do not take into account the magnitude, the 
30 wavelength dependence, or the hemoglobin-composition- 
dependence of the light-scattering effects of unaltered 
whole blood. 



35 



SUMMARY OP THE IKVBHTION 

By making spectrophotometry measurements directly 
on unaltered whole blood and by correcting for the 
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effects of light scattering on the measurements, the 
present invention eliminates the need to hemolyze the 
blood sample and the need to chemically convert the 
various hemoglobin derivatives to a single species. 
5 Thus, the present invention overcomes the previously 
mentioned disadvantages of prior art that employs 
hemolysis only, hemolysis and dilution, or chemical 
conversion. By contrast, the present invention makes 
appropriate measurements of the light scattering by red 

10 blood cells and of other light losses in unaltered, whole 
blood and then uses its assessment of these light losses 
to correct mathematically the measurements of total 
hemoglobin concentration and the concentrations of the 
individual hemoglobin species. By eliminating the need 

15 for a hemolyzing apparatus, the present invention avoids 
the problems caused by hemolysis such as the 
fragmentation of red and white blood cells into light- 
scattering particles. Furthermore, the need for pumps, 
plumbing, and associated control circuitry is eliminated 

20 by use of an inexpensive disposable cuvette. Therefore, 
the problem of clogging of the apparatus is eliminated, 
and the nondestructive optical test performed by the 
present invention preserves the blood sample (with its 
red blood cells intact) for further subsequent analysis. 

25 

The present invention includes the following 
component parts that, employed together and in 
combination, assess the optical transmittance of a sample 
of unaltered whole blood to attain an accurate 
30 measurement of its total hemoglobin concentration, and 

the concentrations of the following pigments of interest: 
bilirubin, oxy-, deoxy-, carboxy-, met-, and 
sulf hemoglobin : 

35 1. A spectrophotometry apparatus (including an 

optical cuvette) in which all optical parameters (for 
example, sample thickness, detector size and shape, 
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sample-to-detector distance, wavelengths , 
and maximum angle of light capture by detector) are 
optimal values so as to minimize the contribution of 
light scattering to the total optical attenuation of 
5 unaltered whole blood and so as to maximize the 
contribution of true optical absorbance. 

2. A method for selecting an optimal set of four 
wavelengths for making the measurements of the 

10 concentrations of oxy-, car boxy-, and methemoglobin with 
the least spectrophotometry error. 

3. A set of seven measuring wavelengths, including 
the four selected by the aforementioned method; wherein 

15 the fifth wavelength is chosen specifically from the 
range of the visible spectrum in which absorption by 
bilirubin is especially large (475 to 500 nm) ; wherein 
the sixth wavelength is chosen specifically from the 
range of the visible spectrum in which absorption by 

20 sulf hemoglobin is especially large (615 to 625 nm) ; and 
wherein the seventh wavelength is chosen from that part 
of the spectrum where absorbance by bilirubin and each of 
the five hemoglobin species is as small as possible in 
comparison with the effects of light scattering, and 

25 small compared with the absorption by at least one of the 
pigments of interest at each of the other six 
wavelengths. The seventh wavelength, the longest 
wavelength, is preferably shorter than 675 nm to avoid 
errors due to the possible presence of indocyanine green 

30 dye in the sample and longer than 660 nm to avoid the 
relatively large absorbance of methemoglobin. 

4. A second method to correct the aforesaid 
hemoglobin and bilirubin concentration measurements for 

35 the effects of light scattering by red blood cells, by 
other specific factors, and by nonspecific. factors in 
unaltered whole blood, wherein the errors due to the 
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t 

various light-scattering losses are treated as functions 
of the wavelength of the measuring radiation and as 
mathematical functions of the concentrations of oxy-, 
car boxy-, met-, deoxy-, and sulf hemoglobin present in 
5 each analyzed sample. 



5. A third method to correct the aforesaid 
concentration measurements for the effects of finite 
bandwidth of the substantially monochromatic measuring 
10 wavelengths, wherein the errors due to finite bandwidth 
of the substantially monochromatic measuring wavelengths 
are treated as functions of the concentrations of oxy-, 
carboxy-, met-, deoxy -r, and sulf hemoglobin present in 
each analyzed sample. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 



20 



Figure 1 is a graph of absorbance spectra for the 
five primary hemoglobin species. 

Figure 2 is a graph comparing the light scattering 
of unaltered whole blood with that of hemolyzed blood. 



Figure 3 is a graph of the theoretical wavelength 
25 dependence of the scattering of light by red blood cells. 

Figure 4 is a graph of the wavelength dependence of 
the light scattering by red blood cells in unaltered 
whole blood in the present invention for varying 
30 hemoglobin compositions of whole blood. 

Figure 5 is a graph of the wavelength dependence of 
the nonspecific light-scattering losses in unaltered 
whole blood for varying hemoglobin compositions of whole 
35 blood. 
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rigure 6 is a block diagram of the present invention 
including the optical geometry of the present invention. 

Figure 7 is a graph of the errors in the 
measurements of relative oxyhemoglobin concentration with 
and without the corrections for nonspecific light 



Figure 8 is a flow chart of the method in accordance 
10 with the present invention to select four principal 
measuring wavelengths. 

Figure 9 is a flow chart of a first method of 
operating the apparatus of Figure 6, in accordance with 
15 the present invention, to correct for light scattering. 

Figure 10 is a flow chart of a second method of 
operating the apparatus of Figure 6, in accordance with 
the present invention, to correct for light scattering. 



Figure lla-e includes a graphical computer program 
listing of an embodiment of the method of Figure 8. 

DBTATliBD DESCRIPTION 

1. Glossary of Terms 



1. Specific Scattering Factor: Known factor 
causing light scattering, and thus light losses, in 
30 unaltered whole blood, the predominant example of which 

is light scattering due to red blood cells (RBCs) . Other 
specific scattering factors include lipid particles, cell 
fragments, microscopic clots, and light-sieving effects 
of sedimented RBCs* 
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2. Nonspecific Scattering Factor: Unknown factor 
causing light scattering, and thus light losses in 
unaltered whole blood 

5 3. Total Light-scattering Losses: The sun of light 

losses due to light scattering by red blood cells, other 
known causes, and unknown causes. In other words, the 
sum of light losses due to Specific and Nonspecific 



10 

4. a. HbC>2 = oxyhemoglobin 

b. HbCO = carboxyhemoglobin 

c. Hi = methemoglobin or hemiglobin 

d. Hb = reduced hemoglobin 
15 e. SHb = sulf hemoglobin 

f . br = bilirubin 

g. Tbr = [br] = total bilirubin concentration 



5. Total Hemoglobin Concentration (THb) : Sum of 
20 the amounts of the various hemoglobin species divided by 
the volume of blood, i.e., the total amount of hemoglobin 
in a mixture of oxy-, carboxy- , met-, deoxy-, and 
sulf hemoglobin divided by the volume of blood. THb = 
[Hb02] + [HbCO] + [Hi] + [Hb] + [SHb]; THb' = [HbC^] + 

25 [HbCO] + [Hi] + [Hb] 



6. Hemoglobin Fraction: Amount of a given species 
(such as oxyhemoglobin) divided by the sum of the amounts 
of oxy-, carboxy- , met-, deoxy-, and sulf hemoglobin , 

30 e.g., fraction of HbC^ = 0.75 

7. Saturation or Relative Concentration: 
Hemoglobin fraction expressed as a percentage, e.g., HfcK>2 
saturation = 75% 

35 

8. Apparent Optical Density (OD) or Apparent 
Absorbance (A) : Logarithm of the total attenuation in 
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light intensity, i.e., OD = A = log 10 (I 0 /I) , where I 0 is 
the incident intensity and I is the intensity transmitted 
through the sample. Note that this is not a re- 
definition of Beer's Law because the apparent optical 
5 density in unaltered whole blood could be due to a 
combination of true Beer's Law absorbance and light 
losses due to scattering. 

9. Error Index: The sum of the absolute values of 
10 the errors that are induced in the measurement of the 
relative concentration of the given hemoglobin species 
(either oxy-, carboxy-, or methemoglobin) due to a 1% 
(relative) change in each of 0Dp... # 0D4* 

15 10. Unaltered Whole Blood: Whole blood that has 

been neither hemolyzed nor diluted. 

2. Optical Geometry* Other Apparatus, and Mode of 
Operation 

20 

Figure 6 shows the presently preferred embodiment of 
the present invention which includes a digital computer 
13; a computer-controlled light source 10 capable of 
sequentially emitting at least seven different 

25 substantially monochromatic wavelengths; an ultra-thin 
disposable optical cuvette 11; a large-area light 
detector 12; and electronic circuitry 14 to amplify the 
detected signal. In operation, a sample of unaltered 
(undiluted and unhemolyzed) whole blood is placed in 

30 cuvette 11 and irradiated with selected wavelengths of 
radiation chosen in accordance with the criteria and 
other principles outlined below. Thus, the normal 
measurement cycle comprises 1) recording the dark current 
from the photodetector when the light source 10 is turned 

35 off, 2) recording the incident light intensity at each of 
the selected wavelengths, 3) placing a sample of 
unaltered (undiluted and unhemolyzed) whole blood in the 
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cuvette 11 and inserting it into the light path, 4) 
recording the light intensity transmitted through the 
sample at each of the selected wavelengths, 5) computing 
the apparent optical densities at each of the selected 
5 wavelengths, 6) correcting the concentration measurements 
for the effects of light scattering by the red blood 
cells and for nonspecific light-scattering losses, 7) 
correcting for the effects of the finite spectral 
bandwidth of the plurality of substantially monochromatic 
10 wavelengths, and 8) reporting the final results to the 
operator. In one embodiment, all measurements are made 
at a substantially constant temperature. 

Following such a measurement cycle, the operator can 
15 remove the disposable cuvette from the instrument and 
discard the sample, or he can withdraw the blood sample 
from the cuvette and subject the unaltered blood sample 
to further hematological analyses. 

20 In one embodiment, light source 10 may be a 

samarium-neon hollow cathode lamp. In another 
embodiment, light source 10 may be a bank of light- 
emitting diodes (LED) arranged so that the emissions of 
each LED pass through a narrow-band optical interference 

25 filter before falling on the sample. In yet another 

embodiment, light source 10 may be a white light source 
such as a halogen lamp optically coupled to a computer- 
controlled monochromator. In the first two embodiments, 
light source 10 may be pulsed at a frequency of 

30 approximately 20 Hz in order to achieve a high signal-to- 
noise ratio by using synchronous detection. 

Disposable cuvette 11 has a short optical path 
chosen to minimize total radiation scattering relative to 
35 absorbance by minimizing the number of scattering 

interactions that a photon is likely to undergo with 
RBCs. Cuvette 11 may have an optical pathlength of 80 to 
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150 ftm. In one embodiment, the optical pathlength is 
90±10 fim. 

Large-area radiation detector 12 is sensitive to the 

t 

5 selected measuring radiation frequencies and may be 
placed within a distance d of from 0 to 10 mm, from 
cuvette 11. Detector 12 has a large active area and is 
placed close to cuvette 11 in order to capture radiation 
scattered at vide angles by the whole blood sample 

10 contained within cuvette 11. In one embodiment, the 

distance d from cuvette 11 to detector 12 is no more than 
1.0 mm, the active area of detector 12 is at least 150 
mm 2 , and preferably at least 600 mm 2 , and the receiving 
half angle a of the detector 12, with respect to the 

15 point where the incident radiation impinges on the 
sample, is at least 30° and preferably at least 70°. 
This optical design minimizes light losses caused by 
light scattered at large angles. Detector 12 may be a 
United Detector Technology detector, PIN-25DP (active 

20 area = 613 mm 2 ) , or may be an SD-1100-11-21-181 detector 
from Silicon Detector Corporation. 



computer 13 may be any type of programmable 
computer, for example, an Apple Macintosh computer. In 
25 order to render the apparatus as portable as possible, 
digital computer 13 may be a programmed microprocessor, 
such as one of the family of Motorola 68000 



30 Thus, the overall strategy of the optical 

configuration of the present invention is to maximize the 
true optical absorbance of unaltered whole blood and to 
minimize the effects of light scattering so that the 
apparent optical density of the sample measured by the 

35 detector is due primarily to absorbance with as small a 
contribution by light scattering as possible. However, 
it should be noted that the optical geometry presented 
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above, when used alone, does not attain commercially 
acceptable accuracy, i.e. it does not measure the 
relative hemoglobin concentrations to an accuracy of 1% 
or less. The optical apparatus of the present invention 
5 achieves commercially acceptable accuracy when used in 
combination with the corrections (described below) that 
correct the hemoglobin concentration measurements for the 
effects of light scattering by the red blood cells, for 
nonspecific light-scattering losses, and for the effects 
10 of the finite spectral bandwidth of the plurality of 
substantially monochromatic wavelengths. 

3. Choice of Principal Four Wavelengths 

15 In the one embodiment of the present invention, six 

of the seven measuring wavelengths are chosen as 
measuring wavelengths for HbC^, HbCO, Hi, Hb, SHb, and 
bilirubin. Four principal wavelengths among the first 
six are chosen to measure HbC^, HbCO, Hi and Hb and are 

20 chosen to minimize the conversion factors by which errors 
in the optical density measurements (ODj, . . • ,OD 4 ) convert 
to errors in the measurements of %HbC>2# %HbCO, and %Hi. 
The method for selecting these four wavelengths is 
described as follows. 

25 

A fixed hypothetical sample of hemolyzed blood whose 
hemoglobin composition is 85% HbC>2 and 15% Hb, whose 

total hemoglobin concentration is 16.11 g/dl (10 mmol/L) , 
is considered in a hypothetical cuvette whose pathlength 
30 is 100 jim. OD£,...,OD4 are defined as the absorbances of 

the hypothetical sample at the four wavelengths \i,...,X 4 , 

thus 

OD l - (0.85 x ejHKE + 0.15 x e^/lO 
OD2 = (0.85 x e 2Hb02 + 0.15 x e^J/lO (1) 
35 OD3 = (0.85 X e 3Hb02 + 0.15 X e 3Hb )/10 

OD 4 = (0.85 x e 4Hb02 + 0.15 x e 4Hb )/10 
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vhere the extinction coefficients shown are the 
millimolar extinction coefficients for the two hemoglobin 
species at the wavelengths Xi,...,\ 4 . Note that once 

\],«..,X 4 are fixed, ZHbOj, %HbCO, %Hi, and %Hb are 

mathematical functions of ODj,,..,OD 4 as shown in the 

following equations: 



ODi/d 


3 lHb02 


e lHbCO 


e lHi 


e lHb 


ODi/d 


: '2Hb02 


e 2HbC0 


e 2Hi 


e Wb 


DD^/d = 


-ZUbOl 


B lHbCO 


e 3Hi 


°3Hb 


0D A /d 


HW>02 


6 4HbCO 


e 4Ht 


e 4Hb 



[HbOj] 

[HbCO] 

[Hi] 

[Bb] 



(2) 



[HbCO] 

[Hi] 

[Hb] 



HHb02 
*2Hb02 
-WbOl 
HHb02 



e lHbCO 
e 2HbCO 
e 3HbCO 
°4HbCO 



e lHi 



e 4JH 



> 














DD 4 /d 



(3) 



THb' = [HbC>2] + [HbCO] + [Hi] + [Hb] ; 
10 %Hb02 = 100 x [Hb0 2 ] / THb'; %HbCO » 100 X [HbCO] / THb' ; 
%Hi = 100 x [Hi] / THb'; %Hb = 100 x [Hb] / THb'; 

Here THb' indicates that SHb is left out of the stun. To 
account for the effects of light scattering when 
15 absorbance is small , Applicants devised the following 
technique: 

Define ODj' « (max {(0.85 x *isb02 + °» 15 x e lHt>) ' 5})/10 
ODj' = (max {(0.85 x ejHbm + 0.15 x eygb) , 5})/l0 
20 OD3' - (max {(0.85 x + 0.15 x e 3Hb ) , 5})/10 

OD 4 ' - (max {(0.85 x e 4HbQ2 + 0.15 x e 4Hb ) , 5})/10 



(The extinction coefficients are the millimolar 
25 absorptivities . ) Thus, whenever the millimolar 

absorptivity falls below 5, then, for the computation of 
ODj', it is set equal to 5. Then for each of the three 
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species (oxy- f car boxy- r and methemoglobin) , an error 
index is computed as the sum of the absolute values of 
the errors that are induced in the measurement of the 
relative concentration of the given hemoglobin species 
5 (either oxy-, carboxy-, or methemoglobin) due to a 1% 
(relative) change in each of ODj' , . . . ,OD 4 ' . 

Thus: 

10 Error Index for oxyhemoglobin = 0.01 

ODj' |d(%Hb02)/d(ODj) |+ 0.01 ODj' | d(%HbC>2) /d(OD2) |+ 0.01 
0D 3 ' |d(%Hb02)/d(OD 3 ) |+0.01 0D 4 ' |d(%Hb02) /d(OD 4 ) |. 

Error Index for carboxyhemoglobin = 0.01 
15 OD 1 '|d(%HbCO)/d(OD 1 ) |+ 0.01 OD2' | d ( %Hbco ) / d ( OD 2> l + °- 01 
003' | d(%HbC0)/d(0D3) 1+0.01 OD 4 ' |d(%HbCO) /d(OD 4 ) | . 

Error Index for methemoglobin = 0.01 
ODj' |d(%Hi)/d(OD 1 ) |+ 0.01 OI^' | d ( %Hi ) / d ( OD 2> l + °- 01 
20 OD3' |d(%Hi)/d(OD3) 1+0.01 0D 4 ' | d(%Hi) /d(0D 4 ) | . 

The derivatives, e.g. d(%HbC0) /d(OI>j) , are taken at 
the fixed values of OD^ r ...,OD 4 determined for the fixed 
hypothetical sample whose hemoglobin is comprised of 85% 
25 Hb02 and 15% Eb. An error criterion for the quadruple of 
wavelengths X 1 ,...,X 4 is then defined as the maximum of 

these three error indices. The mathematical criterion of 
the present invention consists in selecting the unique 
quadruple of wavelengths in the 510-630 nm range which 

30 minimizes the error criterion. The wavelengths are 
searched by 2 nm increments. When the extinction 
coefficient values determined by the inventors are used, 
the final four wavelengths selected by this criteria are 
522, 562, 584 and 600 nm. When the extinction 

35 coefficient values of Zwart et al. are used (A. Zwart, 

E.J. van Kampen and W.G. Zijlstra, Clinical Chemistry 32: 
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972-978, 1986), the final four wavelengths selected by 
this criteria are 518, 562, 580, and 590 nm. 

Referring now to Figure 8, presented is a flow chart 
5 of the method of the present invention to select the 
unique quadruple of wavelengths, in accordance with the 
present invention. 

Beginning in block 81, the Best Error Criterion is 
10 set to 10 10 , and the first quadruple of wavelengths is 

selected. For example, the first quadruple may be 510, 
512, 514, and 516 nm. Control then passes to block 82 
where the Error Index for each of HbC^, HbCO, and Hi are 

computed. Then, control passes to block 83 where the 
15 Error Criterion is set to the absolute value of the Error 
Indices calculated in block 82. 

Then, in decision block 84, the Error Criterion 
calculated in block 83 is compared with the Best Error 

20 Criterion. If the Error Criterion is less than the Best 
Error Criterion, control passes to block 85 where the 
Best Error Criterion is set equal to the Error Criterion 
and the best quadruple of wavelengths is set equal to the 
present quadruple of wavelengths. Control then passes to 

25 block 86. If, on the other hand, decision block 84 

determines that the Error Criterion is not less than the 
Best Error Criterion, control passes directly to decision 
block 86. 

30 In decision block 86, it is determined whether the 

present quadruple of wavelengths is equal to the last 
quadruple of wavelengths, for example, 624, 626, 628 and 
630 nm. If so, control passes to block 88 where the Best 
Error Criterion and best quadruple of wavelengths is 

35 reported. If, on the other hand, block 86 determines 
that the last quadruple has not been processed, control 
passes to block 87 where the next quadruple of 
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wavelengths is formed, and control returns to block 82 
for continued processing* To form the next quadruple, 
only one of the wavelengths at a time is changed by 2 ran. 
This procedure results in the analysis of all possible 
5 combinations of four wavelengths between 510 and 630 nm. 
Any set having duplicate wavelengths is discarded. 

Figures lla-e hereto includes a source code listing 
in LabView graphical language which embodies the method 
10 of the flow chart of Figure 8. In practice, the method 
of Figure 8, either in the form of the source code of 
Figures lla-e, or in an equivalent form, is used to 
program a computer to perform the wavelength quadruple 
selecting method of the present invention. 

15 

4. Corrections for Light Scattering bv RBCs and for 
Nonspecific Light scattering 

As defined in the Glossary, a Scattering Factor is 
20 any specific or nonspecific cause of light scattering, 
i.e. any factor, known or unknown, that causes light 
scattering and hence light losses in unaltered whole 
blood. Examples are red blood cells, lipid particles, 
red blood cell aggregation, clotting, sedimentation, etc. 

25 

According to the present invention, n measuring 
wavelengths are employed to measure k constituent 
components, with n > k, thereby creating an 
overdetermined system of equations with respect to the 

30 chemical compounds being measured. The n * k extra 

equations provide a means by which errors due to n - k 
scattering factors can be compensated. Each specific 
scattering factor in the above definition is assumed to 
contribute in a well-defined way to the apparent 

35 absorbance spectrum of unaltered whole blood placed 

within cuvette C. For a given scattering factor y, the 
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10 



25 



contribution to the apparent optical density of unaltered 
whole blood at wavelength i is denoted OD- iy . 

From these considerations, the following equations 



0D 1 = OD labs + 0D lyl + • • ; + OD lyp 
OD 2 - + ODjyj + . . . + OD^ 



0D n ° 0D nabs + 0D nyl + • • • + OD^ 



where 1, . .., n are wavelength indices and 1, . .., p are 
scattering factor indices. Since there are k constituent 
15 light-absorbing, chemical compounds to be determined, it 
follows that n = p + k. 0D iabs * s the P art of the total 
optical density due to true optical absorption by the k 
light-absorbing compounds, i.e., 

20 0D iahs/ d = e ixl c xl + e ix2 °xZ + ••• + *ixk °xk 

where Xj , • • . , x k denote the k chemical identities of the 
light-absorbing compounds to be measured and c^ , . . . , 
are their concentrations. 



In the context of the present invention, the 
attenuation of light due to the presence of the 
scattering factor y is assumed to obey the relationship 



30 OD 



iy " s iy *y d 

where s y = {s^, S2 y , s^} is the scattering vector 

associated with the scattering factor y, and m^ is the 
magnitude of scattering caused by the scattering factor 
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y. Thus it is possible to write the following matrix 
equation : 



DD2/d 
0D n /d 



e lxl e lx2 
e 2x2 



e L* s lyl s ly2 
e 2xk s 2yl s 2y2 



" s lyp 
• • s 2yp 



3 nxl e nx2 • • • e ruk s nyl s ny2 



• • • S\ 



nyp 



» • • 

c xk 

m yl 
m y2 



yp 



The Applicants have discovered that the scattering 
vectors s y j (where j denotes the scattering factor index) 

5 vary with the hemoglobin species present in the sample 
under test. For example, the scattering vectors s y j are 

different in completely oxygenated blood and blood 
completely saturated with carbon monoxide. Therefore, in 
the present invention, the scattering vectors a^j are 

10 described as functions of the concentrations of the 
hemoglobin species present in the sample under test. 



Since each scattering vector Syj depends on the 
hemoglobin species composition of the sample under 
15 one further distinguishes such scattering vectors as 
SyjfHbC^) , Syj(HbCO) , etc. For an arbitrary sample of 
unaltered whole blood, the appropriate vector is 
computed as 
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Here h is the number of hemoglobin species present among 
the total number k of light-absorbing compounds, and Xj, 
. . . , x h denote the h hemoglobin species being measured . 

5 According to the present invention, the following 

scattering factors have been identified: 1) light 
scattering by red blood cells, a specific scattering 
factor denoted by the subscript RBC, and 2) a nonspecific 
scattering factor, denoted by the subscript NS. The 

10 nonspecific scattering factor, discovered unexpectedly by 
the Applicants, significantly increases the utility of 
the present invention. Data showing the utility of the 
nonspecific scattering factor over and above the 
corrections for light scattering by red blood cells are 

15 given in Figure 7 and Table II following the present 
disclosure of the general method to correct for light 
scattering. 



Disclosed are two embodiments of the corrections for 

20 light scattering: Embodiments A and B. In Embodiment A, 
the following rules hold: only light scattering by the 
red blood cells is considered (so the nonspecific 
scattering correction is ignored) , and the concentrations 
of six compounds are measured, namely HbC^, HbCO, Hi, Hb, 

25 SHb, and br. In Embodiment B, the following rules hold: 
both light scattering by red blood cells and nonspecific 
light scattering are considered, and the concentrations 
of five compounds are measured, namely Hb02, HbCO, Hi, Hb 
and br. Both embodiments use the same set of seven 

30 wavelengths. As explained above, one wavelength (488 nm) 
is chosen at which absorbance by bilirubin is relatively 
large, another wavelength (622 nm) is chosen at which 
absorbance by sulf hemoglobin is relatively large, and 
another wavelength (672 nm) is used because light 

35 scattering by red blood cells is great in comparison with 
absorbance by bilirubin and the hemoglobin species. In 
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addition, the mathematical criterion selected an optimal 
set of four wavelengths (522,562,584, and 600 nm) for 
measuring Hb02, HbCO, Hi with maximum accuracy. If the 
light source is either a battery of light-emitting diodes 

> 

5 combined with optical interference filters or a white 
light source with a monochromator , these seven 
wavelengths can be used in either Embodiment A or B of 

the light-scattering corrections. However, if the 

~> 

samarium-neon hollow-cathode lamp is used, the most 
10 appropriate wavelengths among its emission lines are 

488.4, 520.1, 562.4, 585.2, 597.5, 621.7, 671.7 nm. 

Therefore, this set of wavelengths can be used in both 
Embodiments A and B of the light-scattering corrections. 

15 Once the seven wavelengths have been chosen, values 

for the extinction coefficients for bilirubin and the 
hemoglobin species can be determined. Therefore, the 
following equations obtain in Embodiments A and B of the 
light-scattering corrections . 
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Embodiment A: 

0I W d = e iHb02 1&X>2) + e iHbCO 

+ eiSHb [SHb] + e ibr [br] 



eon, [Hb] 



DD 3 /d 
0D A /d 
DD 5 /d 
0D 6 /d 
ODj/d 



3 lHbQ2 



Hmoi 

9 5Hb02 
3 6Hb02 
*lHb02 



e lHbCO 
e 2HbCO 
e 3HbCO 

G \mco 

e SHbCO 
e 6HbCO 
e lHbCO 



e 2H? 
e 4H! 



B 6Hi 
*11R 



e lHb 
e 2Hb 
e 3Hb 
e 4Hb 
e SHb 
e 6Hb 
e 7Hb 



e lSHb 
B lSHb 
e 3SH& 
e ASHb 
e SSHb 
B 6SHb 
e lSHb 



e lbr 
e 2br 

e 3br 
e 4br 
*5br 
°6br 



s UtBC 
S 2RBC 
S 3RBC 
S 4KBC 
S 5RBC 
SatBC 
SJRBC 



[Sbo 2 i 

[JBbCD] 
[Hi] 
[ffb] 
[SHb] 

[br] 

n KBCscat 



(6A) 



5 THb = [HbQ 2 ] + [HbCO] + [Hi] + [Hb] + [SHb] 
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5 1AS 


5 


S 2NS 


S 


S 3NS 


5 


HNS 


= $ 


*5NS 


> 


*6NS 


3 


S 7NS 


5 



•>4 NS (Hb0 2 ) 
>5Ns(Hb02) 



s 1NS {HbCO) 
S2NslHbCO) 
s 3NS (HbCO) 
s 4NS (HbCO) 
B 5NS {HbCO) 

S 7NS (HbCO) 



s 3NS {Hi) 
s 7NS (Hi) 



s 3NS( Bb 
s 7NS( Hb 



[Hb&i\ /THb 1 

[HbCO] /THb 
[Hi ] /THb 



I 



[Hb]/THb 



(8B) 



wherein [Hbbj] , [HbCO], [Hi], [Hb] , and [SHb] are the 
concentrations of the indicated hemoglobin species in the 
5 sample, ntRBCscat and m NSscat are magnitudes of light 

scattering by the red blood cells and by the nonspecific 
factor, THb = [HbOj] + [HbCO] + [Hi] + [Hb] + [SHb], and 
THb' = [Hb02] + [HbCO] + [Hi] + [Hb] . ODj, . .., ODy are 
the apparent optical densities of a sample of unaltered 

10 whole blood measured at each of the seven wavelengths, 
488.4, 520.1, 562.4, 585.2, 597.5, 621.7, 671.7 nm, 
respectively. The method of finding the correct values 
for the scattering vectors and 8j^ s when the relative 

concentrations (e.g., [Hb02]/THb) are initially unknown 

15 proceeds as follows. Default values for the scattering 
vectors B%g C and S NS are chosen to be the scattering 

vectors 8 RBC (Hb02) and 8 NS (Hb02) , i.e. those values that 
are valid in case the sample is comprised purely of 
oxyhemoglobin . 

20 

With the vectors 8psc(Hb02) and 8 NS (Hb02) placed into 
the matrix equation 6A or 6B above, initial values for 
[HbOj], [HbCO], [Hi], [Hb], [br] , i^RBCscat' ^ C SHb 3 
(Embodiment A) or m NSscat (Embodiment B) are computed. The 
25 initial hemoglobin and bilirubin concentrations are not 
yet accurate enough to report as final results, but they 
are sufficiently close to the true values to compute more 
accurate Srbc and 8 NS vectors for the given blood sample 
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according to equation 7 A or equations 7B and 8B. With 
the newly computed B^ c vector (Embodiment A) , or the 
both the new S^bc and S NS vectors (Embodiment B) placed 
into the matrix equation 6A or 6B above, the quantities 
5 [Hb0 2 ], [HbCO], [Hi], [Hb] , [br], m^c^, and [SHb] or 

m NSscat are recomputed. These quantities are more accurate 
than their initial values because the B^q or both B^q 
and S NS vectors used in obtaining them are more accurate 
than those used in the initial estimates. This procedure 

10 of using the most recently computed hemoglobin 

concentrations to compute more accurate B^qq or both B^ c 
and S NS vectors and in turn using the most recently 
computed Srbc or bot ^ B BBC and S NS vectors to compute more 
accurate hemoglobin and bilirubin concentrations provides 

15 an iterative scheme whereby the hemoglobin and bilirubin 
concentrations quickly converge to their final values. 
In both embodiments A and B, the fourth iteration of the 
computation of these concentrations is sufficiently 
accurate to report. 

20 

Table I shows the values for the scattering vectors 
SrbcO^)' S NS (Hb02), SjiBcOIbCO), S^HbCO) , SreCHb) , and 
S NS (Hb) , which the Applicants have measured and have 
found to be valid for both present Embodiments A and B of 
25 the light-scattering corrections. Figures 4 and 5 show 
the same data graphically at each of the measuring 

wavelengths (488.4, 520.1, 562.4, 585.2, 597.5, 621.7, 
and 671.7 nm) • The present embodiment of the invention 
does not include measured values for the scattering 
30 vectors 8 RBC (Hi) , S NS (Hi) or B mc (SHb) because these less 

common hemoglobin species usually occur in low 
concentrations and because the most common hemoglobin 
species HbO,, HbCO, and Hb usually dominate the blood 
composition in healthy and diseased patients. However, 
35 in Table I are given default values for the vectors 
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SjacCHi), S NS (Hi) and SflBcfSHb) . Moreover, the equations 
are general, so it is within the scope of the present 
invention to include measured values for these scattering 
vectors as well. 

5 

Table I. Scattering Vectors for Various Hemoglobin 
Species 

10 In this table the indices 1 , . . . , 7 refer to the seven 

wavelengths 488.4, 520.1, 562.4, 585.2, 597.5, 621.7, and 
671.7 nm, respectively. 



a. Scattering vectors for HbO^: 

15 S ms (Hb02) 41.33 

(Hb0 2 ) 50 . 00 

SjNsCHbOz) 109.37 

S^sO&Oj) 80.45 

Sjns (HbOj) 3 . 18 

20 SgNsfHbOj) 2.00 

S 7NS (Hb0 2 ) 0.81 

SiKBC( Hbo 2> 3 ' 375 

SajBcflibOj) 2.869 

25 SgnBcdibOj) 2.977 

S 4FBC ( Hbo 2) 3.133 

SfiUBcCHbO^ 1.109 

S 7RBC (Hb0 2 ) 1.036 

30 

b. Scattering vectors for HbCO: 

S 1NS (HbCO) 16.46 

S^sCHbCO) 44.27 

S 3NS (HbCO) 114.30 

35 S 4NS (HbC0) 10.47 
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S SNS (HbCO) 0.90 
S 6NS (HbCO) 1.32 
S^sCHbCO) 1.10 

SiHBcCHbCO) 3.288 

5 S 2RBC (HbCO) 3 . 152 

S 3RBC (HbCO) 3.274 

S 4SBC (HbCO) 2.619 

s 5RBC( Hbc °) 1 ' 618 

10 Syi^cfHbCO) 1.024 

c. Scattering vectors for Hi: 

S 1NS (Hi) 6.216 

SjNsCHi) 6.629 

15 S 3NS (Hi) 7.065 

S 4NS (Hi) 7.200 

S^sCHi) 8.944 

SfiNs(Hi) 21.220 

S ?NS (Hi) 2.800 

20 StKBcCHi) 1.000 

SjRBcCHi) 1.000 

SjpBcCHi) 1.000 

S 4RBc( Hi > 1 - 000 

S 5RBC( Hi ) 1.000 

25 SaacOli) 1.000 

STRBcCHi) 1.000 



vectors for Hb: 

S 1NS (Hb) -38.13 

30 SjNsCHb) -70.24 

S 3NS (Hb) -82.10 

S 4NS (Hb) -81.07 

S^Hb) -41.43 

S 6NS (Hb) -10.81 
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S^sCHb) -4.03 

SiEBC^) 2 - 037 

S 2 RBc( Hb ) 2 - 063 

S 3EBC (Hb) 2.108 

5 S 4KBC (Hb) 2.130 

SfiRBc( Hb ) 1 - 259 

SyRscCHb) 1.092 

10 e. Scattering vector for SHb: 

SiEB C ( SHb > 1 ' 000 

S 2RBc( SHb ) 1.000 

S 3RBC (SHb) 1.000 

S 4RBC (SHb) 1.000 

15 S^cfSHb) 1.000 

S6RB C ( SHb ) 1 - 000 

s 7RBC( SHb ) 1.000 



20 The nonspecific scattering factor (NS) discovered by 

the Applicants provides a refinement of the more coarse 
correction for scattering by the red blood cells (RBC) . 
Table II and Figure 7 show that the nonspecific 
scattering correction improves the method over and above 

25 the accuracy achieved by the RBC correction alone. Table 
II contains the errors in the measurement of unaltered 
. whole blood fully saturated with carbon monoxide with and 
without the nonspecific scattering correction. 
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Table II. Efficacy of the Corrections for Nonspecific 
Light scattering. 

Standard Deviation of %Hb0 2 %HbCO %Hi %Hb 
5 Error 

RBC Correction only 1.41 3.19 2.25 0.50 

RBC & nonspecific 0.67 . 1.28. 0.88 0.53 

correction 

10 

5. Corrections for the Effect of the Finite Spectral 
Bandwidth of the Light Source 

According to the present invention, a method is used 
to correct the measurements of the concentrations of 
Hb0 2 , HbCO, Hi, Hb, SHb, and br for the effects of the 

i 

finite bandwidth of the substantially monochromatic 
measuring wavelengths. This method is not necessary when 
the measuring wavelengths provided by light source B are 
sufficiently monochromatic, as is the case when the 
substantially monochromatic wavelengths are provided by 
isolated emission lines of a hollow cathode lamp. 
However, if the half -intensity bandwidths of the 
measuring wavelengths are greater than 1.0 nm, as they 
would be for a battery of LEDs and optical interference 
filters, this method proves beneficial. 

The corrections for the finite spectral width of the 
light source are necessary because the molar extinction 
30 coefficients of the hemoglobin species vary rapidly with 
wavelength. Thus, when a given spectral width of 
substantially monochromatic light is used, a range of 
(wavelength-dependent) extinction coefficients is 
represented. For fixed spectral bandwidths of the 
35 substantially monochromatic wavelengths provided by light 
' source B and over the range of ODs encountered in 

operation of the apparatus of the present invention, the 
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extinction coefficients of the hemoglobin species are 
linear functions of the concentrations of HbC>2, HbCO, Hi, 
Hb, SHb, and bilirubin present in the sample under test, 
and the corrections in concentrations (given below) are 
5 linear functions of [Hb0 2 ], [HbCO], [Hi], [Hb], [SHb], 

and Tbr. 

When corrections for the finite spectral bandwidth 
of the light source are necessary, the order of operation 

10 of the invention should be first to measure the optical 

densities at each of the seven wavelengths, then to apply 
the corrections for light scattering by the red blood 
cells and for nonspecific light scattering losses, and 
finally to make the corrections for the finite spectral 

15 bandwidth of the light source. 

The method for the compensation of errors due the 
effects of the finite bandwidths of the substantially 
monochromatic wavelengths emitted by light source B 

20 proceeds as follows. After values for the concentrations 
have been corrected for the effects of light scattering 
by the red blood cells and for nonspecific light- 
scattering losses, the corrected concentration values, 
[Hb02], [HbCO], [Hi], [Hb] , [SHb], and Tbr, are then used 

25 to compute spectral bandwidth corrections according to 
the following matrix equations: 
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The final concentrations are then computed according 
to the following formulae: 

[HbO^fa^ = [Hb0 2 ] + [Hb0 2 ] 

correction 

[HbCO]^ = [HbCO] + [HbCO] correcdon 
5 [Hi]^ = [Hi] + [Hi]^^^ 

[Hb] ^1 = [Hb] + [Hb] 

correction 
[SHb]^ = [SHb] + [SHb] 

correction 

[br]^ - [br] + [br] Mrrection 

These matrix equations are completely general in 
that the measurement error that is due to the effects of 
finite bandwidth can depend on the hemoglobin composition 
of the sample under test, mien this method is used, the 
extinction coefficients in matrices 6A and 6B must be the 
extrapolated values for the limit of small 
concentrations. The matrices in equations 9 A and 9B 
depend on the values of the particular measuring 
wavelengths and on the individual bandwidths and spectral 
distributions of the light at each of these wavelengths. 

Referring now to Figure 9, presented is a flow chart 
of the method of operating the apparatus of Figure 6 to 
perform the method of Embodiment A, to correct for light 
scattering in accordance with the present invention. 

Beginning in block 91, the counter n is initialized 
to zero. Then, in block 92, scattering vector Orbc * s 
set equal to Srbc (HbC^) . Then, in block 93 equation 6A 
is used to solve for [HbC^] , [HbCO], [Hi], [Hb], [SHb], 
[br] , and mRBCscar 

Then, control passes to block 94 where the counter n 
is incremented by one. In block 95, it is determined 
whether the counter n is greater than four. If not, 
35 control passes to block 96 where the quantities 

calculated for the concentrations in step 93 are inserted 
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into equation 7A to compute a new value for Srbc 
Control then returns to block 93. If, on the other hand, 
decision block 95 determines that counter n is greater 
than four, control passes to decision block 97 where it 
5 is determined whether the spectral width is greater than 
1 nm. If so, control passes to block 98 where the finite 
bandwidth correction of equation 9A is performed. Then, 
in block 99, the latest values for the concentrations of 
the six constituent components are reported. If, on the 
10 other hand, decision block 97 determines that the 

spectral width is not greater than 1 nm, control passes 
directly to block 99 for report of the concentrations. 

Appendix A includes a source code listing in the C 
programming language which embodies the method of the 
flow chart of Figure 9. In practice, the method of the 
flow chart of Figure 9 embodied in the source code 
program of Appendix A or an equivalent program, is used 
to program computer 13 of Figure 6, to cause the 
apparatus of Figure 6 to perform the method of Embodiment 
A to correct for light scattering. 

* 

Referring now to Figure 10, presented is a method of 
operating the apparatus of Figure 6 , in accordance with 
25 the present invention, to perform the light-scattering 
correction of Embodiment B. 

Beginning in block 101, the counter n is initialized 
to zero. Then, in block 92, scattering vector &rbc * s 
30 set equal to s RBC (Hb0 2 ) and s NS is set equal to s NS (HbC>2) . 
Then, in block 103 equation 6B is used to solve for 
[HbC^], [HbCO] , [Hi], [Hb] , [br] , m^c^, and m NSscar 



15 



20 



35 



Then, control passes to block 104 where the counter 
n is incremented by one. In block 105, it is determined 
whether the counter n is greater than four. If not, 
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control passes to block 106 where the quantities 
calculated for the concentrations in step 103 are 
inserted into equations 7B and 8B to compute a new value 
for B mc . Control then returns to block 103. If, on the 
5 other hand, decision block 105 determines that counter n 
is greater than four, control passes to decision block 
107 where it is determined whether the spectral width is 

■ 

greater than 1 nm. If so, control passes to block 108 
where the finite bandwidth correction of equation 9B is 

10 performed. Then, in block 109, the latest values for the 
concentrations of the six constituent components are 
reported. If, on the other hand, decision block 107 
determines that the spectral width is not greater than 1 
nm, control passes directly to block 109 for report of 

15 the concentrations. 

Appendix B includes a source code listing in the C 
programming language which embodies the method of the 
flow chart of Figure 10. In practice, the method of the 
20 flow chart of Figure 10, embodied either in the source 
code program of Appendix B or an equivalent program, is 
used to program computer 13 of Figure 6, to cause the 
apparatus of Figure 6 to perform the method of Embodiment 
B to correct for light scattering. 

25 

Neither Embodiment A nor B allows the simultaneous 
use of the nonspecific scattering factor and the 
measurement of sulfhemoglobin. A third specific 
embodiment, Embodiment C, which utilizes one more 

30 wavelength than Embodiments A and B, accommodates both of 
these features. Embodiment C uses the same first seven 
wavelengths as do Embodiments A and B. In Embodiment C, 
the following rules hold: both light scattering by red 
blood cells and nonspecific light scattering are 

35 considered, and the concentrations of six compounds are 
measured, namely Hb0 2 , HbCO, Hi, Hb, SHb, and br. 
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In Embodiment C, the eighth wavelength is chosen so 
that the differences e 6Hb02 " e 8Hb02' e 6HbCO " e 8HbCO' e 6Hi ~ 
e 8Hi' e 6Hb " e 8Hb' e 6br " e 8br' e 6RBC " e 8KBO and e 6NS ~ e 8NS' 

are simultaneously small while b 6SS ^ - e 8SHb is quite 

5 large, i.e., in the 635 to 645 nm range. The choice of 
the eighth wavelength, in view of the mathematical 
relationship between a matrix and its inverse, assures 
that the computation of sulfhemoglobin concentration is 
primarily dependent on the 6th and 8th wavelengths and is 
10 affected very little by the other six wavelengths. 
Embodiment C preferably uses 638 nm as its eighth 
wavelength. If the samarium-neon hollow-cathode lamp is 
used, the preferred eighth wavelength among its emission 
lines is 638.3 nm. 
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wherein [HbOj] , [HbCO] , [Hi], [Hb] , and [SHb] are the 
concentrations of the indicated hemoglobin species in the 
sample, niRBCucat and muscat aT& magnitudes of light 
scattering by the red blood cells and by the nonspecific 

5 factor, and THb = [HbO^] + [HbCO] + [Hi] + [Hb] + [SHb] . 

ODj , . . . , ODg are the apparent optical densities of a 

sample of unaltered whole blood measured at each of the 
eight wavelengths, 488.4, 520.1, 562.4, 585.2, 597.5, 
621.7, 671.7, and 638.3 nm, respectively. 

10 

Shown in the following addendum to Table I are the 
additional scattering parameters which are necessary (in 
addition to the parameters in Table I) for use in 
Embodiment C. 

15 

Addendum to Table I . Values for s 8NS ( • ) , s 8RBC ( • ) , and 
s NS (SHb) for Embodiment C. In this addendum the indices 
1,...,8 refer specifically (and only) to the seven 
wavelengths 488.4, 520.1, 562.4, 585.2, 597.5, 621.7, 
20 671.7, and 638.3 nm, respectively. 



s 8NS (Hb02) 
1.379 



s 8NS (HbCO) 
0.682 



s gNS< Hi ) 
22.259 



S 8NS < H3t> ) 
-7.463 



s 8NS( SHb ) 
22.259 



25 s 8RBC( Hb0 2> 
1.051 



SgRBc(HbCO) s 8RBC (Hi) s 8RBC (Hb) s^cfSHb) 
1.117 1.000 1.202 1.000 



s 1NS (SHb) 

S2NS< SHb ) 
30 s 3NS (SHb) 

s 4NS (SHb) 
S 5NS (SHb) 
s 6NS (SHb) 

35 s 8NS (SHb) 



6.216 

6.629 

7.065 

7.200 

8.944 

21.220 

2.800 

22.259 
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6. Compara+ * TO ^ m p lft 

To demonstrate the utility of the present invention 
in comparison with prior art, the Applicants made 
5 measurements of the concentrations of various hemoglobin 
species using a Radiometer OSM3 Hemoximeter and using the 
present invention. The Radiometer OSH3 employs an 
ultrasonic hemolyzer, and is believed to embody the 
optical geometry of Johansen et al. (U.S. Patent 
10 3,972,614), and the method and apparatus of Lundsgaard 
(U.S. Patent 4,997,769) to correct for the turbidity of 
hemolyzed blood. 

In normal operation, the OSM3 first ultrasonically 
15 hemolyses each sample, then makes measurements of the 

sample's optical absorbance at multiple wavelengths, and 
finally uses the method of Lundsgaard to correct for the 
residual turbidity of the hemolyzed sample. To determine 
whether the method and apparatus of Lundsgaard could 
20 accurately measure hemoglobin concentrations in unaltered 
whole blood, the Applicants analyzed blood samples with 
the OSM3 operating normally and then analyzed the same 
blood samples with the 0SM3's ultrasonic hemolyzer 
disabled • 

25 

Table III shows that the Radiometer OSM3, believed 
to use the method and apparatus of Lundsgaard (U.S. 
Patent 4,997,769), gives completely erroneous results 
when the blood sample is not hemolyzed. For example, 

30 measurements of the total hemoglobin concentrations with 
and without hemolysis differ by as much as 7.4 grams/dl. 
Similarly, it is physically impossible for the relative 
concentration of oxyhemoglobin to be 400% or for the 
concentration of carboxyhemoglobin to be negative (-99%). 

35 Thus, it is readily apparent that the method and 
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apparatus of Lundsgaard fail to yield even plausible 
results on unaltered, whole blood. 
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To evaluate the present invention's ability to make 
accurate measurements directly on unaltered, whole blood, 
the Applicants employed the herein described methods and 
apparatus that constitute the present invention and made 
5 measurements of the concentrations of various hemoglobin 
species in blood samples before and after hemolysis. The 
results appear in Table IV. 
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Table IV shows that the measured total hemoglobin 
concentrations with and without hemolysis typically 
differ by only a few tenths of a gram/dl. Similarly, the 
measurements of the relative concentrations of oxy-, 
5 car boxy-, and methemoglobin are consistently in close 

agreement in all of the samples, in summary, the results 
shown in Tables III and IV demonstrate conclusively that 
the present invention enables one skilled in the art to 
make accurate measurements of the hemoglobin 
10 concentrations directly in unaltered, whole blood; 

whereas the method and apparatus of Lundsgaard fail to 
obtain satisfactory results. 

While the present invention has been described in 
15 connection with preferred embodiments, it will be 

understood that these embodiments are only examples, and 
should not be considered a limitation of the present 
invention. 
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APPENDIX A 

COPYRIGHT, 1992, BOARD OF REGENTS, THE UNIVERSITY OF TEXAS SYSTEM 
/*********************^ 

/* FUNCTIONS: 1) void slope_comp_LED0 */ 
I* 2) void OD_to_7sat_LED0 */ 

/* 3) void slope_recur_LED() */ 

^*********************&^ 

include "G\A^^VOX^C_Pll(X^TYPE^AVOX_TYP.C ,, 
#include "C\ATG^AVOX>C_PROG\DEIAAVOX_DEF.C" 
#include ' ' Q\ATG\A VOXSC_PROG\MATH\MATH_DEF. C ' ' 
include "C\ATGV^VOX\C_PROGP«DR\AVOXJHDR-C" 

extern const vector Insats_vec; 
extern const matrix Identity_mtx; 
extern const matrix Scat_mix; 
extern const matrix Sulf_mtx; 
extern const matrix THbslp_mtx; 

i(XinsaLXscat,Xsulf^CTHb) 



{ 



vector *Xinsat; I* Xinsat (Input) = Initial %Sats* */ 
vector *Xscau /*Xscat (Output) = 'Scat Slopes (Vector)' */ 
vector *Xsulf; /*Xsulf (Output) = "Sulf Slopes (Vector)' */ 
vector *XTHb; /* XTHb (Output) = 'THb Slopes (Vector)' */ 



AtimesX((matrix *)Scat_mtx,XinsaUCscat); 
AtimesX((matrix *)Sulf_mu^Gnsat^sulf); 
AtimesX((matrix *)THbslp_mtx,Xinsat^THb); 
return; 

} f* end slope_comp_LED0 */ 

voidOD_to_7sat_LED(AExtC2,ACalS^\Csats^CODs^CSHbs^Csats^ScatCF^satslP, 

xTotHb^cTHbc3nSlp,EnTHb,EnVer) 
matrix *AExtC2; /* AExtC2 (Input) = "Extinction Coeffs.2' */ 
matrix *ACalS; I* ACalS (Input) = 'Calibration Sats' */ 
matrix *ACsats; /* ACsats (Input) = 'Computed Sats* */ 
vector *XODs; I* XODs (Input) = 'ODs' */ 
vector *XSHbs; I* XSHbs (Output) = "SHb Correction Factors' */ 
vector *Xsats; /*Xsats (Output) = 'SATS' */ 
vector *XScatCF; I* XScatCF (Output) = 'Scattering Corr'n Fac' */ 
vector *XsatslP; I* SsatslP (Output) = 'SATS.lst-Pass' */ 
ljeal *xTotHb; /*xTotHb (Output) = Total Hb' */ 
l_realxTHbc; /*xTHbc (Input) = THb constant" */ 
boolean EnSlp; /*EnSlp (Input) = 'Enable Slopes' */ 
boolean EnTHb; /*EnTHb (Input) = "Enable THb Slopes' */ 
boolean EnVer, /*EnVer (Input) = "Black-Old version; Wht-'*/ 

{ 

matrix Atmp 1 AtmpZ^AtmpS; 

vector Xtmp 1 ,Xtmp2^Ctmp3^tmp4^tmp5 r Xtmp6^Ctmp7 ,Xtmp8 ; 
l_real xtaipl^tmp2,xtmp3^tmp4^tmp5,xtrnp6,xto^ 
Lreal sulf_thb_cutoff, sulf_slp_cutoff; 
boolean enable_slp, enable_thb_slp; 

Ainverse((matrix *)A(^,(matrix *)Atmpl); 
AtimesB((matrix *)Atmpl, (matrix *)ACsats,(mattix *)Atmp2); 
AtimesB ((matrix *)Atmp2,(matrix *)AExtC2,(matrix *)Atmpl); 
Atranspose((matrix *)Atmpl, (matrix *)Atmp2); 
Ainverse((matrix *)Atmp2,(matrix *)Atmp3); 
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/* switch and calculate Y from XODs - do not change Atmp3 yet */ 
xtmpl = getXvalue((vector *)XODs,0); 
xtmp2 = x_abs(xtmpl); 
xtmp3 = xtmpl / xtmp2; 

Xtimesx((vector *)XODspctmp3,(vector *)Xtmpl); /* Xtmpl = Y */ 

AtimesX((matrix *)Atmp3,(vector *)Xtmpl,(vector *)Xtmp2); /* may discard Atmp3 now */ 

Xsubset((vector *)Xtmp2,0 f 5,(vector *)Xtmp3); 

xsum = Xsum((vector *)Xtmp3); 

Xdivx((vector *)Xtmp2 > xsum,(vector *)Xtmp 1); 

eep_set_var((uJong)&enable_slp, ,, Enable_Slq)es n ); 

eep_set_var((uJong)A^nable_thb_slp, M EnableJIHb_Sl^^ 

sulf_thb_cutoff = 0_real)-L0; 
sulf_slp_cutoff = (Lreal)-l-0; 

/* slope_recur_LED #1 */ 
I* parameter = variable 10 */ 
/* Xsatsf = Xtmpl */ 
/*Xconc =Xtmp2 */ 
/* Xscat = Xtmp3 (output-discarded) */ 
/* Xsulf = Xtmp4 (output-discarded) */ 
/* Xsats = XtmpS */ 
/* xsulflHb = sulf_thb_cutoff */ 
/* xsulfslp = sulf_slp_cutoff */ 
/*xTHbc = xTHbc */ 
/* EnSlp = enable„slp */ 
I* EnTHb = enable_thb_slp */ 

slope_recur_LED((vector *)Xtmpl, 

(vector *)Xtmp2, 
(vector *)Xtmp3, 
(vector *)Xtmp4, 
(vector *)Xtmp5, 
(l_real)sulf_thb_cutoff, 
Q_ieal)sulf_slp_cutoff, 
(Lreal)xIHbc, 
(boolean)enable_slp, 
(boolean)enable_thb_slp); 



Xtimesx((vector *)Xtmp5,(l_real)100A(vector *)XsatslP); /* define SATS.lsLPass */ 

/* slope_recur_LED #2 */ 
/* parameter = variable 10 */ 
/* Xsatsf =Xtmp5 */ 
/♦Xconc =Xtmp2 */ 
/♦Xscat =Xtmp7 */ 
/♦Xsulf =Xtmp8 */ 
/♦Xsats = Xtmpl */ 
/* xsulflHb = sulfjhb_cutoff*/ 
/* xsulfslp = sulf_slp_cutoff*/ 
/* xTHbc = xTHbc */ 
/* EnSlp = enable„slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recur_LED((vector ♦)Xtmp5, 

(vector ♦)Xtmp2, 
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(vector *)Xtmp7, 
(vector *)Xtmp8, 
(vector ♦)Xtmpl, 
(l_real)sulf_thb_cutoff, 
Q_real)sulf_slp_cutoff, 
G_real)xTHbc, 
(boolean)enable_slp> 
(boolean)enable_thb_slp); 

/* slope_recur JJED #3 */ 
/* parameter = variable 10 */ 
/*Xsatsf =Xtmpl */ 
/*Xconc =Xtmp2 */ 
/* Xscat = Xtmp3 (output-discarded) */ 
/* Xsulf = Xtmp4 (output-discarded) */ 
/* Xsats = XtmpS */ 
/* xsulfTHb = sulf_thb_cutoff ♦/ 
/* xsulfslp = sulf_slp_cutoff */ 
/* xTHbc = xTHbc */ 
/* EnSIp = enable_slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recur_LED((vector *)Xtmpl, 

(vector *)Xtmp2, 
(vector *)Xtmp3, 
(vector *)Xtmp4, 
(vector *)Xtmp5, 
Q_real)sulf_thb_cutoff, 
O_real)sulf_slp_cutoff, 
(LreaOxTHbc, 
(boolean)enable_slp, 
(boolean)enable_thb_slp); 

/* slope_jecur_LED #4 */ 
/* parameter = variable IO */ 
/*Xsatsf =Xtmp5 */ 
/♦Xconc =Xtmp2 */ 
/♦Xscat =Xtmp7 */ 
/♦Xsulf = Xtmp8 */ 
/♦Xsats = Xtmpl */ 
/♦ xsulflHb = sulf_thb_cutoffV 
/* xsulfslp = sulf_slp_cutoff*/ 
/♦xTHbc =xlHbc ♦/ 
/* EnSIp = enable_slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recur_LED((vector ♦)Xtmp5, 

(vector *)Xtmp2, 
(vector *)Xtmp7, 
(vector *)Xtmp8, 
(vector *)Xtmpl, 
(l_real)sulf_thb_cutoff, 
(l_real)sulf_slp_cutoff, 
(LreaDxTHbc, 
(boolean)enable_slp, 
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(boolean)enable_thb_slp); 



Xsubset((vector *)Xtmp5,0,4,(vector *)Xtmp3); 

Xtimesx((vector *)Xtmpl,0_real)100.0,(vector *)Xsats); /* define Xsats */ 

Xclear((vector *)Xtmp6); 

Xtmp6[0] = (Lrcal)0.020; 

Xtmp6[l] = (Lreal)0.020; 

Xtmp6[2] = (Lreal)O.OOO; 

Xtmp6[3] = (Lreal)aOl 6; 

xtmp4 = XYJnpioduct((vector *)Xtmp3,(vector *)Xtmp6); 
xtmpl =Xtmp2[6]; 
xtmp7 = xtmpl * xtmp4; 
xtmp8 = Xtmp2[4]; 
xtmp9 = xtmp7 + xtmp8; 

xtmp4 = xtmp9 / xsum; /* xtmp4 gets compared to -1.0 */ 
Xclear((vector *)Xtrap7); 
Xtmp7[0] = 0_real)-7.50; 
Xtmp7[l] = 0_real)-2.50; 

= (Lreal) 0.00; 
= (Lreal) 12.00; 
xtmp5 = XY_inproduct((vector *)Xtmp3,(vector *)Xtmp7); 
Xclear((vectar *)Xtmp8); 
Xtmp8[0] = 0-ieal)-0.280; 
Xtmp8[l] = (Lreal)-0.280; 
Xtmp8[2] = (Lreal) 0.000; 
Xtmp8[3] = (LrealMUSO; 

xtmp6 = XY_inproduct((vector *)Xtmp3,(vector *)Xtmp8); 
if(EnSlp = FALSE) { 

xtmpS = (l_real)0.0; 

xtmp6 = (l_real)0.0; 
} /* end if Enable Slopes is FALSE */ 



Xtmp7[2 
Xtmp7[3 



if(xtmp4 > (l_real>L0){ xtmpS = (Ijreal)O.O; } 

xtmp4 = xtmpS * xtmp9; 

xtmp2 = xtmp6* xtmpl; 

xtmp3 = xtmp4 + xtmp2; 

xtmp4 = xsum +xtmp3; 

*xTotHb = xTHbc * xtmp4; 

return; 

} f* end ODjto_7saL.LED0 */ 



void slope_recurJLED(Xsatsf, Xconc, Xscat, Xsulf, Xsats, xsulfTHb, 

xsulfslp, xTHbc, EnSlp, EnTHb) 
vector *Xsatsf; f* Xsatsf (Input) = 'Sats (fractional)' */ 
vector *Xconc; /*Xconc (Input) = 'Concentrations 1 */ 
vector *Xscat; /* Xscat (Output) » 'Scat slopes* */ 
vector *Xsulf; /*Xsulf (Output) = *Sulf slopes 1 */ 
vector *Xsats; /* Xsats (Output) = 'SATS' */ 
Lreal xsulflHb; /* xsulfTHb (Input) = 'SulfTHb Cut-off */ 
Lreal xsulfslp; /* xsulfslp (Input) = 'Sulf Slope Cut-off */ 
Lreal xTHbc; /* xTHbc (Input) = THb constant 1 */ 
' boolean EnSlp; /*EnSlp (Input) = 'Enable Slopes' */ 
boolean EnTHb; /* EnTHb (Input) = 'Enable THb Slopes' */ 

{ 

vector Xtmpl ,Xtmp2 t Xtrnp3,Xtirip4; 
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Ureal xtmpl t xtmp2,xtmp3,xsum; 
long intl; 

/* function call: slope_compJLEDO */ 
/* parameter = variable I/O */ 
/* Xinsat = Xsatsf I */ 
/* Xscat = Xscat O */ 
/* Xsulf = Xsulf O */ 
/* XTHb = Xtmpl O */ 

slope_comp_LED((vector *)Xsatsf,(vector *)Xscat,(vector *)Xsulf,(vector *)Xtmpl); 

if(EnlHb = TRUE) [ xtmpl = (Lreal)LO; } 
else{ xtmpl = (Lxeal)0.0; } 

Xtimesx((vector *)Xtmpl f xtmpl .(vector *)Xtmp2); 

Xsubset((vector *)Xconc,0,5,(vector *)Xtmp3); 

xsum = Xsum((vector *)Xtmp3); 

xtmpl = xTHbc * xsum; 

xtmp2 = xtmpl - (Lreal)13.5; 

Xtimesx((vector *)Xtrap2,xtmp2,(vector *)Xtmp 1); 

Xdivx((vector *)Xtmp 1 ,0_ieal) 1 00.0,(vector *)Xtmp3); /* save Xtmp3 */ 

if(EnSlp = TRUE){ xtmpl = (Lreal)l.O; } 
else{ xtmpl = Q_real)0.0; } 

Xtimesx((vector *)Xscat,xtmpl, (vector *)Xtmpl); 
xtmp2 = getXvalue((vector *)Xconc,6); 
Xtimesx((vector *)Xtmpl f xtmp2, (vector *)Xtmp2); 
XplusY((vector *&tmp2,(vector *)Xconc,(vector *)Xtmpl); 
xtmp3 = getXvalue((vector *)Xtmpl,4); 
xtmp2 = xtmp3 / xsum; 
if(xtmp2 > xsulfTHb){ intl = 5; } 
else{ intl =4; } 

if(xtmp2>xsulfslp){ xtmp2 = (Lreal)O.O; } 
else{ xtmp2 = (Lreal)l.O; J 

Xtimesx((vector *)Xsulf pctmp 1 /vector *)Ximp4); 
Xtimesx((vector *)Xtmp4pctmp2,(vector *)Xtmp2); 
Xtimesx((vector *)Xtmp2ixtmp3,(vector *)Xtmp4); 
XpIusY((vector *^tmpl,(vector *)Xtmp4,(vector *)Xtmp2); 
Xsubset((vector *)Xtmp2,04ntl t (vector *)Xtmpl); 
xtmp3 = Xsum((vector *)Xtmpl); 
Xdivx((vector *)Xtmp2^tmp3,(vector *)Xtmp4); 
xtmpl = Xtmp2[5]; 
xtmp2 = xTHbc* xtmpl; 

Xtmp4[5] = xtmp2; /* place value into vector */ 
XplusY((vector *)Xtmp3,(vector *)Xtmp4,(vector *)Xsats); 
return; 

} /* end slope_recur_LEDO */ 



/* END FILE: MTX_COMP.C */ 
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APPENDIX B 

COPYRIGHT, 1992, BOARD OF REGENTS, THE UNIVERSITY OF TEXAS SYSTEM 
y**************** ********** 

/* FUNCTIONS: 1) void slope_comp_LED() */ 

I* 2) void OD_to_7saLLED0 */ 

I* 3) void slope_recur_LED () */ 

include T:\ATG\AVOXSC.PROG\TYEE^VOX_.TYP.C M 
include "C:\ATG^VOX\C_PROG^DEF\AVOX_DEF.C" 
#include "C:\ATG\A VOXNC_PROG\MAra\MATFLPEF.C" 
include "C:\ATG\AVOX\C^PROGSHDR\A VOX_.HDR.C n 

extern const vector Insats_vec; 
extern const matrix Identity_mtx; 
extern const matrix Scat_mtx; 
extern const matrix Sulfjntx; 
extern const matrix THbslp_mtx; 

void slope^comp^LEDfXinsauXscauXsulf^XTHb) 

vector *Xinsat; /* Xinsat (Input) = Initial %Sats' */ 
vector *Xscau /* Xscat (Output) = 'Scat Slopes (Vector)' */ 
vector *Xsulf; /*Xsulf (Output) = 'Sulf Slopes (Vector)' */ 
vector *XTHb; /* XTHb (Output) = THb Slopes (Vector)' */ 

{ 

AtimesX((matrix *)Scat_mtx T Xinsat > Xscat) ; 
AtimesX((matrix *)Sulf_mtx,XinsauXsulf); 
AtimesX((matrix *)THbslp_mtx > XinsauXTHb); 
return; 

} /* end slope_comp_LED0 */ 

void OD_toJ7saU^(AExtQ^CalS AC 

xTotHbptTHbc3nSlpJEnTHb,EnVer) 
matrix *AExtC2; /* AExtC2 (Input) = 'Extinction Coeffs.2' */ 
matrix *ACalS; /*ACalS (Input) = 'Calibration Sats' */ 
matrix *ACsats; /*ACsats (Input) = 'Computed Sats' */ 
vector *XODs; /*XODs (Input) ='ODs' */ 
vector *XSHbs; /*XSHbs (Output) = 'SHb Correction Factors' */ 
vector *Xsats; /*Xsats (Output) = 'SATS' */ 
vector *XScatCF; /* XScatCF (Output) = 'Scattering Corr'n Fac/ */- 
vector *XsatslP; /* SsatslP (Output) = 'SATS.lstPass* */ 
Lieal *xTotHb; /*xTotHb (Output) = Total Hb' */ 
Lreal xTHbc; /*xTHbc (Input) = THb constant' */ 
boolean EnSlp; /* EnSlp (Input) = TEnable Slopes' */ 
boolean EnTHb; /*EnTHb (Input) = 'Enable THb Slopes' */ 
boolean EnVer, /*EnVer (Input) = 'Black-Old version; Wht-*/ 

{ 

matrix Atmpl ,Atmp2Atmp3; 

vector Xtmpl JCtmp2 t Xtmp3^Xtmp4^Ctmp5 > Xtmp6 ,Xtmp7,Xtmp8; 
Lreal xtmp 1 ,xtmp2,xtmp3 ,xtmp4 ,xtmp5 ,xtmp 6»xtmp7 »xtmp 8 »xtmp9 ,xsum; 
Lreal sulf_thb_cutoff, sulf_sIp_cutoff; 
boolean enable_slp, enable_thb_slp; 

Ainverse((inatrix *)ACalS, (matrix *)Atmpl); 
AtimesB((matrix *)Atmpl, (matrix *)ACsats t (matrix *)Atmp2); 
AtimesB((matrix *)Atmp2,(matrix *)AExtC2,(matrix *)Atmpl); 
Atranspose((matrix *)Atmpl, (matrix *)Atmp2); 
Ainverse((matrix *)Atmp2,(matrix *)Atmp3); 
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/* switch and calculate Y from XODs - do not change Atmp3 yet */ 
xtmp 1 = getXvalue((vector *)XODs,0); 
xtmp2 = x_abs(xtmpl); 
xtmp3 = xtmpl / xtmp2; 

Xtimesx((vector *)XODs,xtmp3,(vector *)Xtmpl); /* Xtmpl = Y */ 

AtimesX((matrix *)Atmp3,(vector *)Xtmpl,(vector *)Xtmp2); /* may discard Atmp3 now */ 

Xsubset((vector *)Xtmp2,0,5,(vector *)Xtmp3); 

xsum = Xsum((vector *)Xtmp3); 

Xdivx((vector *)Xtmp2,xsum,(vector *)Xtmpl); 

eep„set_va^((uJong)&enable - slp t ,, Enable_Slopes ,, ); 

eep_set„var((uJong)toi^^^ 

sulf_thb_cutoff = (l_real)1.25; 
sulf_slp_cutoff= (l_real)1.25; 

/* slope_recurJLED #1 */ 
/* parameter = variable IO */ 
/* Xsatsf = Xtmpl */ 
/♦Xconc =Xtmp2 */ 
/* Xscat = Xtmp3 (output-discarded) */ 
/* Xsulf = Xtmp4 (output-discarded) */ 
/*Xsats =Xtmp5 */ 
/* xsulflHb = sulf_thb_cutoff */ 
/* xsulf sip = sulf_slp_cutoff */ 
/*xTHbc =xTHbc */ 
/* EnSlp = enable_slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recurJJED((vector *)Xtmpl, 

(vector *)Xtmp2 t 
(vector *)Xtmp3, 
(vector *)Xtmp4, 
(vector *)Xtmp5, 
O_real)sulfjhb_cutoff, 
O_real)sulf_slp_cutoff, 
(Lreal)xTHbc, 
(boolean)enable_slp, 
(boolean)enable_tbb_slp); 

Xtimesx((vector *)Xtmp5,(l^_real)100.0,(vector *)XsatslP); /* define SATS.lsLPass */ 

/* slope_recurJJED #2 */ 
/* parameter = variable IO */ 
/* Xsatsf =Xtmp5 */ 
/* Xconc = Xtmp2 */ 
/* Xscat =Xtmp7 */ 
/* Xsulf =Xtmp8 */ 
/*Xsats = Xtmpl */ 
/* xsulfTHb = sulf_thb_cutoff*/ 
/* xsulfslp = sulf_slp_cutofF*/ 
/*xTHbc =xTHbc */ 
/* EnSlp = enable_slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recur_LED((vector *)Xtmp5, 

(vector *)Xtmp2, 
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(vector *)Xtmp7, 
(vector *)Xtmp8, 
(vector *)Xtmpl, 
O_ieal)sulf_thb_cutoff, 
(l_real)sulf_slp_cutoff, 
(l_real)xTHbc, 
(boolean)enable_slp, 
(boolean)enable_thb_slp); 

/* slope^.recur_LED #3 */ 
/* parameter = variable IO */ 
I* Xsatsf =Xtmpl */ 
/* Xconc = Xtmp2 */ 
/* Xscat = Xtmp3 (output-discarded) */ 
/* Xsulf = Xtmp4 (output-discarded) */ 
/*Xsats =Xtmp5 ♦/ 
/* xsulfTHb = sulf_thb_cutofif */ 
/* xsulfslp = sulf_slp__cutoff */ 
/* xTHbc = xTHbc */ 
/* EnSlp = enable_slp */ 
/* EnTHb = enable_thb_slp */ 

slope_recurJUED((vector *)Xtmpl, 

(vector *)Xtmp2, 
(vector *)Xtmp3, 
(vector *)Xtmp4, 
(vector *)Xtmp5, 
O_ieal)sulf_thb_cutoff, 
O_ieal)sulf_slp_cutoff, 
(Lreal)xlHbc, 
(boolean)enable_slp, 
(boolean)enable_thb__slp); 

/* slope_recur JJED #4 */ 
/* parameter = variable IO */ 
/* Xsatsf =Xtmp5 */ 
/♦Xconc =Xtmp2 */ 
/♦Xscat =Xtmp7 */ 
/♦Xsulf = Xtmp8 */ 
/♦Xsats =Xmq>l */ 
/* xsulflHb = sulLAb.cutoff*/ 
/♦ xsulfslp = sulf_^)_cutoff*/ 
/♦ xTHbc =xTHbc ♦/ 
/* EnSlp = enable_slp ♦/ 
/* EnTHb = enable_thb_slp ♦/ 

slope_recurJLED((vector ♦)Xtmp5, 

(vector ♦)Xtmp2, 
(vector *)Xtmp7, 
(vector *)Xtmp8, 
(vector ^Xtmpl, 
(Ueal)sulf_thb_cutoff, 
O_jeal)sulf_slp_cutoff, 
(Lreal)xTHbc, 
(boolean)enable_slp, 
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(boolean)enable_thb_slp); 
Xsubset((vector *)Xtmp5,0,4,(vector *)Xtmp3); 

Xtimesx((vector *)Xtmpl,0_real)100.0,(vector *)Xsats); /* define Xsats */ 

Xclear((vector *)Xtmp6); 

Xtmp6[0] = 0_ieal)0.020; 

Xtmp6[l] = (l_real)0.020; 

Xtmp6[2] = (l_real)0.000; 

Xtmp6P] = (l_real)0.016; 

xtmp4 = XY_inproduct((vector *)Xtmp3,(vector *)Xtmp6); 
xtmp 1 = Xtmp2[6]; 
xtmp7 = xtmpl * xtmp4; 
xtmp 8 = Xtmp2[4]; 
xtmp9 - xtmp7 + xtmp8; 

xtmp4 = xtmp9 / xsum; /* xtmp4 gets compared to 1.25 */ 

Xclear((vector *)Xtmp7); 

Xtmp7[0] = 0_real)-7.50; 

Xtmp7[l] = (l_real)-2.50; 

Xtmp7[2] = (ljreal) 0.00; 

Xtmp7[3] = (l_real)12.00; 

xtmp5 = XYJnpioduct((vector *)Xtmp3,(vector *)Xtmp7); 
Xclear((vector *)Xtmp8); 
Xtmp8[0] = 0 jeal)-0.280; 
Xtmp8[l] = (i_real)-0-280; 
Xtmp8[2] = (Lreal) 0.001 
Xtmp8[3] = 0_Jeal)-0.150; 

xtmp 6 = XY_inproduct((vector *)Xtmp3,(vector *)Xtmp8); 
if (EnSlp = FALSE) { 

xtmp5 = (Lreal)O.O; 

xtmp6 = (Lreal)0.0; 
} /* end if Enable Slopes is FALSE */ 

if(xtmp4 > (Lreal) 1 .25) { xtmp5 = (l_real)0.0; ) 
xtmp4 = xtmpS * xtmp9; 
xtmp2 = xtmp6 * xtmpl; 
xtmp3 = xtmp4 + xtmp2; 
xtmp4 = xsum + xtmp3; 
*xTotHb = xTHbc * xtmp4; 
return; 

} /* end OD_to_7sat_LED0 */ 



I slope_recur_ J LED(Xsatsf, Xconc, Xscat, Xsulf, Xsats, xsulfTHb, 

xsnlfslp, xTHbc, EnSlp, EnTHb) 
vector *Xsatsf; /* Xsatsf (Input) = 'Sats (fractional) , */ 
vector *Xconc; /* Xconc (Input) = ■Concentrations' */ 
vector *Xscat; /* Xscat (Output) = 'Scat slopes' */ 
vector *Xsulf; /* Xsulf (Output) = 'Sulf slopes' */ 
vector *Xsats; I* Xsats (Output) = 'SATS' */ 
Lreal xsulfTHb; I* xsulfTHb (Input) = 'Sulf THb Cut-off */ 
Lreal xsulfslp; /* xsulfslp (Input) = 'Sulf Slope Cut-off */ 
Lreal xTHbc; /* xTHbc (Input) = THb constant" */ 
boolean EnSlp; I* EnSlp (Input) = "Enable Slopes' */ 
boolean EnTHb; /* EnTHb (Input) = "Enable THb Slopes' */ 

vector Xtmpl,Xtmp2,Xtmp3,Xtmp4; 
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Lreal xtmpl,xtmp2,xtnip3,xsum; 
long intl; 

I* function call: slope_compJ-ED0 */ 
/* parameter = variable I/O */ 
I* Xinsat = Xsatsf I */ 
I* Xscat = Xscat O */ 
I* Xsulf =Xsulf O */ 
/*XTHb = Xtmpl O */ 

slope_comp_LED((vector *)Xsatsf,(vector *)Xscat,(vector *)Xsulf,(vector *)Xtmpl); 

if(EnTHb = TRUE) { xtmpl = (l_real)1.0; } 
else{ xtmpl - OUreaI)0.0; } 

Xtimesx((vector *)Xtmpl,xtmpl,(vector *)Xtmp2); 

Xsubset((vector *)Xconc,0,5,(vector *)Xtmp3); 

xsum = Xsum((vector *)Xtmp3); 

xtmpl = xTHbc * xsum; 

xtmp2 = xtmpl - fl_ieal)13.5; 

Xtimesx((vector *)Xtmp2,xtmp2,(vector *)Xtmpl); 

Xdivx((vector *)Xtmp l,(Lreal)l 00 A(vector *)Xtmp3); /* save Xtmp3 */ 

if(EnSlp = TRUE){ xtmpl = 0_real)L0; } 

else{ xtmpl = (l_real)0X); } 

Xtimcsx((vector *)Xscat t xtmp 1 ,(vector *)Xtmpl); 

xtmp2 = getXvalue((vector *)Xconc,6); 

Xtimesx((vector *)Xtmp 1 > xtmp2,(vector *)Xtmp2); 

XplusY((vector *)Xtmp2,(vector *)Xconc,(vector *)Xtmpl); 

xtmp3 = getXvalue((vector *)Xtmpl,4); 

xtmp2 =s xtmp3 / xsum; 

if(xtmp2 > xsulfTHb){ intl = 5; } 

else{ intl =4;} 

if(xtmp2 > xsulf sip) { xtmp2 = (l_real)0.0; } 

elsc{ xtmp2 = (Lreal)l.O; } 

Xtimesx((vector *)Xsulf ,xtmp 1 ,(vector *)Xtmp4); 
Xtimesx((vector *)Xtmp4^tmp2,(vector *)Xtmp2); 
Xtimesx((vector *)Xtmp2ixtmp3,(vector *)Xtmp4); 
XplusY((vector *)Xtmpl,(vector *)Xtmp4 t (vector *)Xtmp2); 
Xsubset((vector *)Xtmp2,04ntl, (vector *)Xtmpl); 
xtmp3 = Xsum((vector *)Xtmpl); 
Xdivx((vector *)Xtmp2 f xtmp3,(vector *)Xtmp4); 
xtmpl = Xtmp2[5]; 
xtmp2 = xTHbc * xtmpl; 

Xtmp4[5] = xtmp2; /* place value into vector */ 
XplusY((vector *)Xtmp3,(vector *)Xtmp4,(vector *)Xsats); 
return; 

} /* end slope_recur_LEDO */ 



/* END FILE: MTX.COMP.C */ 
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CLAIM8: 



1. A method of determining the concentrations of a 
5 plurality of constituent components of unaltered whole 
blood, including: 

generating a plurality of substantially 

monochromatic radiation wavelengths, each 

10 wavelength of an absorbance subset of said 

plurality of wavelengths having been selected 
by their ability to distinguish the constituent 
components and having been selected to minimize 
the effects of radiation scattering and to 

15 maximize radiation absorbance by said 

constituent components, and each wavelength of 
a scattering subset of said plurality of 
wavelengths having been selected to maximize 
the effects of radiation scattering by 

20 unaltered whole blood relative to the effects 

of radiation absorbance by unaltered whole 
blood ; 

irradiating a sample of unaltered whole blood with 
25 said plurality of radiation wavelengths, 

through a depth of said sample chosen to 
minimize radiation scattering by unaltered 
whole undiluted blood; 

30 detecting intensities of said radiation wavelengths, 

after passing through said depth of said 
sample, at a distance from said sample, and 
over a detecting area, both chosen to minimize 
the effects of radiation scattering by 

35 unaltered whole blood on the determination of 
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concentration of said constituent components; 
and 

calculating concentrations of said plurality of 
5 constituent components of said sample of 

unaltered whole blood, based upon detected 
intensities of each of said plurality of 
radiation wavelengths, and based upon 
predetermined molar extinction coefficients for 
10 each of said constituent components at each of 

said radiation wavelengths of said absorbance 
subset . 

15 2. The method of claim 1, wherein said depth of said 
sample is in the range of 80 to 150 micrometers. 

3. The method of claim 2, wherein said depth of said 
20 sample is approximately 90 micrometers. 

4. The method of claim 1, wherein said detecting area 
is at least approximately 150 square millimeters. 

25 

5. The method of claim 4, wherein said detecting area 
is at least approximately 600 square millimeters. 

30 

6. The method of claim 1, wherein said distance from 
said sample is within the range of 0 to 10 millimeters. 

35 7. The method of claim 6, wherein said distance from 
said sample is approximately 1 millimeter. 
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8. The method of claim 1, wherein said step of 
detecting is performed over a half-aperture angle of 
radiation emanating from said sample of at least 
approximately 30 degrees. 

5 

9. The method of claim 8 f wherein said step of 
detecting is performed over a half-aperture angle of 
radiation emanating from said sample of at least 

10 approximately 70 degrees. 

10. The method of claim 1, further comprising: 

15 correcting said calculated concentrations of 

constituent components for the effects of 
finite spectral bandwidth of the substantially 
monochromatic wavelengths on the extinction 
coefficients of each constituent component. 

20 

11. The method of claim 1, said plurality of constituent 
components including HbO^, HbCO, Hi and Hb, said method 
further comprising: 

25 

before said generating step, selecting four 

radiation wavelengths by computing an error 
index for each of HbC^, HbCO and Hi as the sum 
of the absolute values of the errors that are 
30 induced in the measurement of relative 

concentrations of Hb0 2 , HbCO and Hi due to a 
change in optical density measurements; and 



35 



selecting a quadruple of radiation wavelengths 
having minimum error indices. 
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12. The method of claim 11, said quadruple of radiation 
wavelengths each being within the range of 510 to 630 
nanometers . 

5 

13. The method of claim 12, said quadruple of radiation 
wavelengths comprising 522, 562, 584 and 600 nanometers. 

10 14. The method of claim 12, said quadruple of radiation 
wavelengths comprising 518, 562, 580 and 590 nanometers. 

15. The method of claim 12, said quadruple of radiation 
15 wavelengths comprising 520.1, 562.4, 585.2 and 597.5 

nanometers . 

16. The method of claim 12, said constituent components 
20 further including bilirubin, said method further 

comprising: 

before said generating step, selecting a radiation 
wavelength within the range of 475 to 500 
25 nanometers as the radiation wavelength for the 

measurement of bilirubin. 

17. The method of claim 16, said radiation wavelength 
30 for the measurement of bilirubin being 488.4 nanometers. 

18. The method of claim 12, said constituent components 
further including sulf hemoglobin, said method further 

35 comprising: 
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before said generating step, selecting a radiation 
wavelength within the range of 615 to 625 
nanometers as the radiation wavelength for the 
measurement of sulf hemoglobin. 



19. The method of claim 18, said radiation wavelength 
for the measurement of sulf hemoglobin being 621.7 
nanometers . 



20. The method of claim 1, further comprising: 

correcting said calculated concentrations of 
15 constituent components for the effects of light 

scattering by red blood cells. 



21. The method of claim 20, said correcting step 
20 comprising, correcting said calculated concentrations of 
constituent components as a function of the relative 
concentrations of the constituent components. 



25 22. The method of claim 21, said correcting step further 
comprising: 

iteratively determining a red blood cell scattering 
vector for the particular composition of the 
30 whole blood sample being analyzed; and 

using said red blood cell scattering vector to 

correct said calculated constituent component 
concentrations • 



35 
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23. The method of claim 1, further comprising: 

correcting said calculated constituent component 

concentrations for the effects of non-specific 
5 light scattering. 



24. The method of claim 23, said correcting step 
comprising, correcting said calculated concentrations of 
10 constituent components as a function of the relative 
concentrations of the constituent components under 
consideration • 



15 25. The method of claim 24, said correcting step further 
comprising: 

iteratively determining a non-specific scattering 
vector for the particular composition of the 
20 whole blood sample being analyzed; and 

using said non-specific scattering vector to correct 
said calculated constituent component 
concentrations • 

25 

26. The method of claim 1, further comprising, 
correcting said calculated concentrations of constituent 
components for the effects of light scattering by red 
30 blood cells and for the effects of non-specific light 



27. The method of claim 26, said correcting step 
35 comprising, correcting said calculated concentrations of 
constituent components as a function of the relative 
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concentrations of the constituent components under 
consideration. 

5 28. The method of claim 27, said correcting step further 
comprising: 

iteratively determining a red blood cell scattering 
vector for the particular composition of the 
10 whole blood sample being analyzed; 



iteratively determining a non-specific scattering 
vector for the particular composition of the 
whole blood sample being analyzed; and 

15 

using said non-specific scattering vector and said 
red blood cell scattering vector to correct 
said calculated constituent component 
concentrations . 



29. The method of claim 28, said plurality of 
constituent components including HbC^, HbCO, Hi and Hb, 
said method further comprising: 

25 

before said generating step, selecting four 

radiation wavelengths by computing an error 
index for each of HbC^, HbCO and Hi as the sum 
of the absolute values of the errors that are 
30 induced in the measurement of relative 

concentrations of Hb0 2 , HbCO and Hi due to a 
change in optical density measurements; and 

selecting a quadruple of radiation wavelengths 
35 having minimum error indices. 
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30. The method of claim 29, said quadruple of radiation 
wavelengths each being within the range of 510 to 630 
nanometers • 



5 

31. The method of claim 30, said constituent components 
further including bilirubin, said method further 
comprising: 

10 before said generating step, selecting a radiation 

wavelength within the range of 475 to 500 
nanometers as the radiation wavelength for the 
measurement of bilirubin. 

15 

32. The method of claim 31, said constituent components 
further including sulf hemoglobin , said method further 
comprising: 

20 before said generating step, selecting a radiation 

wavelength within the range of 615 to 625 
nanometers as the radiation wavelength for the 
measurement of sulf hemoglobin. 



33. The method of claim 32, further comprising, before 
said generating step, selecting a radiation wavelength 
within the range of 635 to 645 nanometers as an 
additional radiation wavelength for the measurement of 
3 0 sulfhemoglobin • 
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